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Abstract 



In this article we suggest a fast multi-scale edge-detection scheme for medical ultrasound 

^>- . signals. The edge-detector is based on well-known properties of the continuous wavelet trans- 

r J I form. To achieve both good localization of edges and detect only significant edges, we study the 

• . maxima-lines of the wavelet transform. One can obtain the maxima-lines between two scales 

f_\ ' by computing the wavelet transform at several intermediate scales. To reduce computational 

effort and time we suggest a time-scale filtering procedure which uses only few scales to connect 

modulus-maxima across time-scale plane. The design of this procedure is based on a study of 

^ . maxima-lines corresponding to edges typical for medical ultrasound signals. This study allows 

\^ ' us to construct an algorithm for medical ultrasound signals which meets the demand for speed, 

OO . but not on expense of reliability. 

The edge-detection algorithm has been applied to a large class of medical ultrasound sig- 
nals including tumour-, liver- and artery-images. Our results show that the proposed algorithm 
g : effecfvely detect. „,ai» reat.tes .„ »ch signals, including edges with low conttast. 

1 Introduction 

^ . Recent advances vi^ithin medical research calls for new techniques to improve the quality of med- 

H I ical ultrasound signals. Ultrasound is versatile and non-invasive compared to many other medi- 

cal visualization-systems. Unfortunately ultrasound signals are often heavily contaminated vi^ith 
speckle-noise. In addition, edges in such signals often have low contrast. This makes medical ultra- 
sound signals difficult to interpret and analyze. Improving the quality of such signals is therefore 
desirable. 

In this article we study a fast wavelet-based edge-detector scheme for medical ultrasound sig- 
nals. This scheme is based on properties of the continuous wavelet transform. We refer to [1] and 
[2, Ch. 6] for an introduction to the wavelet transform as well as its application in edge-detection. 
The continuous wavelet transform can be used to study local structures - such as edges - at different 
scales in a signal. In particular, if the wavelet is given by the 1st derivative (1-D) or gradient (2-D) 
of the Gaussian-function, the local extrema (modulus-maxima) of the wavelet transform at a scale 
correspond to edge-points in the signal. 

In order to explain the idea behind our investigation we will give a brief survey of the now clas- 
sic Canny edge-detector [3] . The Canny edge-detector computes the wavelet transform of a signal 
at a single scale. If the modulus of the wavelet transform at a modulus-maximum is larger than 
a threshold-value its spatial location represents an edge. At a fine scale the Canny edge-detector 
provides fine details, but is sensitive to noise. At a coarse scale one can detect low-contrast edges, 
but the positional accuracy of a modulus-maximum is poor compared to the location of the edge. 
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At each scale the wavelet transform yields a different description of the edges in the signal. It 
is clear that for a noisy signal - such as a medical ultrasound signal - there is no single scale which 
provides a categorically correct representation of the edges. To deal with this ambiguity we - as 
many others e.g. [4, 5, 6, 7, 8, 9] - suggest to use the wavelet transform at several scales. The idea 
is to study modulus-maxima at coarse scales to detect the edges, and use a fine scale to localize 
them. 

The main purpose of this article is to address two issues which frequently occur when using 
multi-scale wavelet transform to detect edges. These are; 

1. Which modulus-maxima at different scales correspond to the same edge in the underlying 
signal, and 

2. how can one effectively decide which modulus-maxima correspond to the same edge. 

As the scale decreases each modulus-maximum belongs to a curve (maxima-line) in the time-scale 
plane. As s— s-O each maxima-line propagates towards an edge in the signal. The maxima-lines 
provide a natural identification of which modulus-maxima corresponds to the same edge in the 
signal. 

To decide which modulus-maxima belong to the same maxima-line one can compute the wavelet 
transform at several close scales and gradually trace the modulus-maxima towards fine scales [4, 5]. 
However, the wavelet transform computed at two close scales carries almost the same information, 
and do not contribute with new information with regards to the edge-detector scheme studied in 
this article. Computing the wavelet transform at several scales will therefore yield redundancy, as 
well as making the procedure computationally expensive. 

Our aim is to decide which modulus-maxima belong to the same maxima-line by computing the 
wavelet transform at only a few 'different' scales. To achieve this it is essential to understand the 
behavior of the maxima-lines in the time-scale plane. This depends on several factors, such as the 
intensity and mutual position of the edges. In general the ambiguity introduced by using multiple 
scales is inescapable. The goal is not to eliminate the ambiguity, but rather manage it and reduce 
it where possible. 

For this purpose we study typical maxima-lines corresponding to medical ultrasound signals. It 
happens that local intensity-changes in medical ultrasound signals often can be approximated by 
a set of model-edges. In this article we refer to these as patterns. In particular, the maxima-lines 
of each pattern have their special behavior in the time-scale plane. 

To the best of our knowledge the idea of using model-edges to study properties of the wavelet 
transform in the time-scale plane is due to [10] who created an 'elementary' set of model-edges. 
Since then this idea has been used by a number of researchers e.g. [11] and more recently [12]. 

We propose a new time-scale filtering procedure for medical ultrasound signals. This procedure 
accommodates the time-scale behavior of both the maxima-lines and the wavelet transform corre- 
sponding to the set of patterns. The procedure uses a decision-function which analyzes a decay- 
and a distance-criterion to decide which modulus-maxima belong to the same maxima-line. The 
procedure is constructed to decide this by computing the wavelet transform at only a sparse set of 
scales. 

The set of patterns serves as test signals for the time-scale filtering procedure. By analyzing 
the accuracy of the decision-function on the maxima-lines corresponding to each pattern we study 
the reliability of the procedure. We complement this with a statistical study on medical ultrasound 
signals. 

We also study an edge-detector scheme for medical ultrasound images which is based on the 
time-scale filtering procedure. The edge-detector uses a perceptual criteria similar to [6] to find the 



significant edges in an image. This method has shown promising results for detecting and locahzing 
(low-contrast) edges in medical signals. 

Several have suggested edge-detectors based on the multi-scale wavelet transform. We refer 
to [13] and references therein for a general exposition of edge-detectors. In [14, 15, 12] edges are 
characterized by estimating the decay of the wavelet transform along each maxima-line. In [6, 9] 
they use multi-scale wavelet coefficients to detect and localize significant edges, and in [16, 17] the 
authors suggest a procedure using automatic scale selection. Statistical methods have been used in 
[18], and in [7] a method based on the multi-scale Edge-flow vector field have been proposed. 

The article is organized as follows. In the next section we provide a brief reminder of the wavelet 
transform and its use to detect edges. We also discuss general properties of maxima- lines related to 
edge-detection. In Sect. 3 we propose and study a time-scale filtering procedure for medical ultra- 
sound signals. We conclude Sect. 3 with a discussion of a space-scale filtering procedure for images. 
This construction is based on our analysis of the 1-D time-scale filtering procedure. In Sect .4 we 
present the edge-detection algorithm. In the final part of the note we present experimental results 
of both time-scale filtering procedure and edge-detection algorithm. These results are based on a 
study of both phantom and actual medical ultrasound images. 

2 Wavelet transform. 

2.1 Wavelet transform. 

We consider the following well-known smoothing- function and wavelet; 

The wavelet transform is defined as e.g. in [2]; 

Wf{u,s)= f f(t)^^C-^^)dt. 

The candidate edges of / correspond to the points (u, s) in the time-scale plane where |VF/(-, s)\ has 
a local maximum. These points are referred to as modulus-maxima, in what follows abbreviated 
mod-max. The set of mod-max of Wf{u, s) at a scale s > will be denoted Aif{s); 

M.f{s) — {m : uisalocal— maximumof |VK/(-, s)|}. 

For 2-D signals (images) J. Canny introduced an edge-detector in [3]. This was later formulated 
in wavelet terminology by e.g. [15]. The smoothing function and (directional) wavelets are given 

by; 

e{x,y) =e{x)e{y), 

V{x,y) ='ip{x)e{y), ipy{x,y) = e{x)ip{y), 

where 9{-) and ^(•) are defined as above. Given a function / in two variables, let fs{x,y) = 
s~^f{s~^x,s~^y). The 2-D wavelet transform is defined as; 

The vector Wf is the direction of maximal change of the smoothed image f*0s- Candidate edges 
are the points [{u,v),s^ in the space-scale plane where |M^/((n, f), s) | has a local maximum in 



p. NES 



the direction of maximal change. The set of mod-max at scale s will in what follows be denoted 
Mfis); 

A4f{s) = \^{u,v) : (u,t;) isadirectionallocalmaximumof M^/((-, •), s)}. 

2.2 Maxima-lines. 

With respect to the wavelet considered in this article it is well-known that every mod-max of A^/(s) 
belongs to some curve 

i = {{u,s) ■.u = iais), a£A}, 

where A is some index-set. These curves are called maxima- lines. For wavelets equal a derivative 
of the Gaussian-function every maxima-line ia{s) propagates towards finer scales [2, page 178], see 
fig.l for an example. In particular, if no maxima-lines propagate towards a point in an interval it 
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Figure 1: (a) A 1-D ultrasound signal, (b) The maxima-lines of Wf{u,s). 

contains no edges [2, Th. 6.5]. 

It will be assumed that a mod-max at a scale corresponds to only one mod-max at any finer 
scales. In theory a maxima-line may bifurcate, but the probability of bifurcating in real signals is 
negligible. 

To illustrate this we consider the function f(t)=p(t)+Ap{t—l) where A>0 and p(t) denotes the 
Heaviside-function (we remind the definition of p{t) in Sect. 3. 2). The mod-max of Wf{u,s) are 
found by solving 
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There exists a s* > (see (3)) such that this equation has one and only one solution for every 
s > s* and three solutions (two modulus-maxima and one local minimum) for s < s* . It is clear 
that the solution curves (maxima-lines) bifurcate if and only if In^ = i.e. A = 1. For all other 
values of A the maxima-lines do not bifurcate. Our study of maxima-lines in medical ultrasound 
signals indicates that bifurcations are rare in such signals. 

We will also assume that two mod-max on the same maxima-line correspond to the same 
underlying phenomenon in the signal, and that the accuracy of the localization of the mod-max 
increases as s decreases. These assumptions corresponds to the classical identity- and localization- 
assumption of Witkin [5] . 

3 Space-scale filtering. 

The main purpose of this article is to study how one can connect mod-max between two (perhaps 
rather distinct) scales. Let ia, d ^ Ahe the maxima- lines of Wf{u, s), and 



Mfis^) = {iais^)■.a£A}, i = l,2 



be the niod-max at two scales. If S2 > si the number of elements in Mf{s2) is less or equal 
the number of elements in A4f(si) since new maxima-lines may appear between S2 and si. The 
question is how one can decide which mod-max of A4/(si) belong to the same maxima-line as a 
mod-max of Aif{s2), and which belong to a maxima- line which appears between S2 and si. 

At fine scales one can decide which mod-max belong to the same maxima-line by considering 
their spatial distance. At coarse scales one can consider the logarithmic-decay of the wavelet 
transform. In addition the sign of the wavelet transform along a maxima-line can not change. 

In the next section we suggest a time-scale filtering procedure for medical ultrasound signals. 
This procedure is based on a decision- function P(-, •) which analyzes both a distance- and a decay- 
criterion to decide which mod-max belong to the same maxima-line. We design the function to 
favor decay at coarse scales and distance at fine scales. In Sect. 3. 2 we introduce a set of model- 
edges, called patterns, which serves as test-signals for the decision-function. These patterns model 
local intensity-changes in medical ultrasound signals. In Sect. 3.3 we study the reliability of the 
time-scale filtering procedure with respect to each pattern and their respective maxima- lines. 

3.1 Decision- function and time-scale filtering procedure. 

Let S2, si be two scales S2 > si, and let (n, S2) and (m, si) be two mod-max. We define the decision- 
function corresponding to scales S2 and si as; 

P(n,m) = A{n,m)D{n,m)Sign{n,m), (1) 

where; 

A(n,m) = exp 

r)C„ ™\ - exn ( - In l^'^^"'''^^! ]n~^ £2 _ l 
u{n,m) — (ix.p\ m 1^^^^ ,,^^1 m ^^ ^ 




Sign(n,m)=<; " '^ ^'^'^\Wfim,s^)) 
otherwise. 

The essence of the procedure is to connect mod-max which maximize the decision- function P(-, •). 
For instance, if {n,S2) and {m,si) both belong to the same maxima-line ia say, then we want 
P(ra, m) > P(n, u) for all u G Mf{si), u ^ m. The control-parameter a in the expressions of A and 
D controls the mutual weighting between the distance- (A) and decay-criterion (D). A negative a 
favors the decay-criterion, while a positive favors distance. By studying the procedure for various 
choices of a we find a weighting between distance and decay suitable for our type of signals. 

In the following sections we will analyze the decision function P(-,-) and the corresponding 
time-scale filtering procedure. This analysis is based on a set of model-edges which model local 
intensity changes in ultrasound signals. 

3.2 Patterns 

To investigate whether the suggested procedure successfully connects mod-max between two scales 
we have studied a collection of model-edges, called patterns. These patterns - six in total - model 
typical transitions between tissues in medical ultrasound signals, fig. 2. 
Let px denote the translated Heaviside-function; 

, . ( I t>x 
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Figure 2: Illustration of how one can locallly model edges in medical ultrasound signals with the 
patterns from Sect. 3. 2. 



The patterns can be modeled as; 



Pattern 1 
Pattern 2 
Pattern 3 
Pattern 4 
Pattern 5 
Pattern 6 



f{t)=pa{t)~Bpp{t), 

f{t)^Po{t)+Ap,{t), 

f{t)^Po{t)+Ap,{t)-Bpp{t), 

f{t)^Po{t)+Api{t)-Bp0{t), 

f{t)=po{t)-Bpp{t)+Api{t), 



B,f3>0, 

A>1, 

A,p>l, B*>B>0, 

P>1,B>0, A*>A>0, 

1>I3>0,A,B>0. 



The six patterns each model a particular type of edge. Pattern 1 models an isolated step-edge, 
that is, if the intensity at the boundary of the object changes at a single jump. Pattern 2 models a 
thin object or an impulse, while Pattern 3 models a boundary where the intensity changes at two 
isolated jumps. Pattern 4 and pattern 5 model an object with boundary similar to Pattern 3, but 
where the intensity decreases shortly after the initial intensity-changes. This occurs for instance if 
there is a 'void' within the object. If the object has a 'hole' near the boundary one can model this 
with Pattern 6. 

Pattern 4 and pattern 5 differ by the time-scale behavior of the maxima-lines, fig.3(h,j). Given 
A, f3 (or B, (3) the point of transformation between the patterns are given by B* (or A*). B* or A* 
can be found by solving; 

dWf{u,s) 
du 







d^Wfiu,s) 

d^Wf{u,s) 
d-ip 




0. 



(2) 



Fig. 3 shows typical examples of the patterns and their respective maximal-lines. Every pattern 
can after a suitable normalization and translation be assumed to have an edge at with intensity 1. 
Pattern 3-6 have (after rescaling) an additional edge at 1 with (positive) intensity A, and pattern 
2,4-6 have an edge of (negative) intensity B at /3 ^ 0,1. The six patterns have one (pattern 1), two 
(pattern 2,3) and three (pattern 4-6) maxima-lines, in what follows denoted Iq, (.\ and ^/j. 

In fig. 3 the six patterns and characteristic examples of their maxima-lines are displayed. As 
indicated in fig.3h, one can expect that distance between mod-max is a weak criterion for pattern 
4. Rather surprisingly this is also the case for pattern 3. Observe that the maxima-lines of these 
'simple' patterns show a complexity similar to those in medical ultrasound signals. 

A common factor of pattern 3-6 is that there exists a scale where the maxima-line £o (pattern 
3,4,6), li (pattern 5) or ij^ (pattern 6) first appears. This scale is found by solving the equations; 



dWf(u,s) 

d^Wf(u,s) 
du'^ 




0. 



The unique s > for which (3) has a solution is denoted s* 
filtering procedure is often most critical near s* . 



(3) 



As it happens, the suggested time-scale 



(a) Pattern 1: p{t) 



(b) Maxima-line of pattern 1 



(c) Pattern 2: p{t)-Bpp(t) 





(d) Maxima-lines of pattern 2 



(e) Patterns: p{t) + Api{t) 
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(f ) Maxima-lines of pattern 3 



(g) Pattern 4: p{t)+Api{t)~Bpp{t) 
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(h) Maxima-lines of pattern 4 



(i) Pattern 5: p(t)+Api(t)-Bp^(t) 
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(j) Maxima-lines of pattern 5 
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(k) Pattern 6: p[t)-Bpfi{t)^Apx(t) 
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(1) Maxima-lines of pattern 6 



Figure 3: The dictionary of patterns (model edges) used in our analysis, (a) step-edge (c) impulse- 
edge (e) staircase-edge (g)-(i)-(k) various triplet edges 



3.3 Evaluation of the time-scale filtering procedure. 

The purpose of the forthcoming analysis is to study for which patterns (i.e. for which A,B, 13) the 
time-scale filtering procedure suggested in Sect. 3.1 works. We also study the procedure for different 
control-parameters a. 

Let as before Iq, l\ and l^ denote the maxima-lines in the patterns from Sect. 3. 2 and fix two 
scales S2 and s\. We chose in this study 81 = 82/2^ i.e. the scales are dyadic distributed. In what 
follows let P(-, •) be the decision-function corresponding to 82 and si defined in Sect. 3.1. 

It is straightforward to see that the procedure works for the first two patterns, hence we focus 
our analysis on patterns 3—6. The decision-function P(-, •) needs to fulfill one inequality for each 
maxima-line. Namely; the decision-function connects two mod-max ^1(52) and ^i(si) from the 
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(a) 



(b) 



Figure 4: (a) The graph of Qi{A) and (b) the graph of Qo{A) for 1 < ^ < 1.5. The 
cyan, yehow, blue, red, green curves correspond to P(n,m) defined by (1) with respectively 
a = {—1/2,0,1/2,1,3/2} in the expression of P. The dashed line corresponds to P = D and 
dash-dotted line corresponds to P = A. 



maxima-line ii if; 

P{h(s2)Ji{si))-P{h{s2),ej{si)) >0 fori = 0,/3. 

Similarly the mod-max io{s2) and io{si) are connected if; 

P(4(s2),4(si))-P(4(s2),^j(si)) >0 forj = l,/?. 



(4) 



(5) 



Note that for pattern 3 the inequalities (4) and (5) have to hold only for j = and j = 1 respectively 
since ifs does not exist. Finally, for pattern 4-6 the suggested procedure connects two mod-max 
from the maxima-line Ir if; 



P{^fs{s2)j0{si)) -P{ef,(s2)Jj{si)) > forj-0. 



(6) 



Inequality (6) always holds for pattern 4-6 since Sign(^^(s2),^j(si)) =0 only if j ^ (3 (since 
Wf[ii3{s),s) <0 and Wf[ij{s),s) >0, j = 0, l). We will therefore focus the analysis on the in- 
equalities (4) and (5). 

For pattern 3,4,6 the decision-function successfully connects any two mod-max if both inequal- 
ities (4) and (5) hold for all scales such that £o exists. Consequently (4) has to hold for every 
0<si<s* and (5) for all 0<S2<s*. For pattern 5 both (4) and (5) have to hold for all scales such 
that ii is well-defined. Thus (4) has to hold for every 0<S2<s* and (5) for all 0<,si<s*. 
Comment: The preceding analysis of the decision-function P(-, •) focuses on patterns with values 
of A,B,f3 typical for our specific application (locally model medical ultrasound signals). Values of 
A,B and/or j3 near the 'transformation-limit' (value of A,B,f3 where a pattern transforms into 
another) will not be studied. In practice one can often consider such 'critical' configurations by 
another pattern. 

Pattern 3: First we consider the decision-function and corresponding time-scale filtering 

procedure with respect to pattern 3. Numerical analysis suggests that the time-scale filtering pro- 
cedure is most likely to fail at the largest scale where both ii and £o exist. We can therefore simplify 
the analysis by considering the left-hand-side (Ihs.) of (4) with si = s*, S2 = 2si and Ihs. of (5) 
with S2 = s* , si = S2/2. Note that the Ihs. of (4) and (5) - denoted Qi{A) and Qo{A) - are now 
functions depending only on A. In particular, if both Qi{A) > and QoiA) > the suggested 
procedure holds. 

Fig.4 displays the graphs of respectively Qi{A) and Qo{A) for P(-, •) given by (1) with control- 
parameters a={— 1/2, 0, 1/2, 1,3/2} (cyan, yellow, blue, red, green). The dashed and dash-dotted 
lines correspond to the 'limit cases' resp. P(n,m)=D{n,m) and P(n,m)=A{n,m). 

From fig. 4a one can see that along the maxima-line £1 the decay-criterion (dotted line) is op- 
timal for connecting mod-max; if P(n,ni) = D(n,m) then Qi(A) > (i.e. (4) holds) for all A > 1. 
Along Iq one need to have A> 1.3, fig. 4b. On the contrary; if one considers the distance-criterion 



(dash-dotted line) it suffices to have A>1 at fine scales, but needs to have ^ > 1.3 at coarse scales. 
Consequently one must have A> 1.3 in order to connect mod-niax between any two scales S2 = 2si 
and < si < s* for pattern 3 if either P(n, m) = D(n, m) or P(n, m) = A{n, m). 

Combining distance and decay increase the ^-interval for which the time-scale filtering proce- 
dure works, fig. 4. If e.g. P(n, m) is given by (1) with q = (yellow line) then it suffices to have 
A > 1.1 to connect mod-max between (any) two scales in pattern 3. If a = 1/2 (blue line) then both 
Qi(A)>0 and Qq{A)>^ for ah ^> 1.16. 

Pattern ^; For pattern 4 there are three maxima-lines, io, ii and ip. As illustrated in 




.B 2 n ^i, 2.6 2.8 




(a) Q = 1/2 



(b) a = -1/2 



Figure 5: Numerical analysis of Q\ for pattern 4. The figures display level curves of the surface 
(5i(y4,S,/3) = for P(n,m) given by (1) with (a) a = 1/2 and (b) a = -1/2. NOTE that in (a) 
we consider 1.1 < /3 < 3 and in (b) 1.1 < /3 < 2. 



(a) a = -1/2 

Figure 6: Numerical analysis of Qo for pattern 4. The figure displays some level-curves of the 
surface Qo{A,B, 13) = for P(n, ?tt,) given by (1) with a = —1/2. 

fig.3h the edge at /3 'repels' the maxima-lines corresponding to the edges at and 1. The result is 
that the distance between io and ii gets smaller, hence weakening the distance-criterion for pattern 
4. 

For pattern 4 the time-scale filtering procedure successfully connect mod-max across scales if 
both (4) and (5) holds. For this purpose let Qi{A, B, (5) denote the minimum of the Ihs. of (4) with 
respect to S2 = 2si, < si < s*, j = and < B < B*, where B* is given by (2) (recall that pattern 
4 transforms into pattern 5 for B > B*). Similarly let Qq{A, B ^ f3) denote the minimum of the Ihs. 
of (5) with si = 2s2,0<S2<s*, j = l andO<5<S*. li both QiiA,B,f3)>0 and Qo{A,B,l3)>0 
the procedure successfully connects any two mod-max with respect to pattern 4. 

The additional variables B and /3 make the analysis of Qi and Qq more complicated for pattern 
4-6 than for pattern 3. Rather than considering the graphs of Qi and Qq we study the (level-curves 
of the) surfaces Qi{A, B, 13) = 0, i = 0,1. Numerical analysis indicates that both Qi and Qq increase 
if: (i): the jump at location 1 increases, (ii): the jump at (3 decreases, or (Hi): the distance from 
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1 to /3 increases. In particular, if Qi{Aq, Bq, j3o) = then Qi{A, B, P) >0 for all A > Aq, Bq> B >0 
and /? > /3o, and similarly for Qq. 

As with pattern 3; different values of the control-parameter a 'interpolates' P(-, •) given by (1) in 
between the 'limit' cases P(n, m) = D(n, m) and P(n, m) = A(n, m). Fig. 5 displays the level-curves 
of the surface Qi{A,B, (3) = with respect to P(n,m) given by (1) with a = it 1/2 (note that 
the y-axis differs in the figures). Fig. 6 displays the level-curves of Qo{A,B, f3) = for a = —1/2. 
Numerical studies indicate that Qo{A, B, I3)>0 if a = 1/2, except possibly near critical values of 
A,B,j3 such as e.g. /3~ 1 (which are uninteresting to study for our purpose). 

As expected; increasing the influence of the decay-criterion increases the accuracy of the time- 
scale filtering procedure for pattern 4. The following table illustrates this for a few concrete values 
of A and /3. The table shows the maximal B-interval such that the time-scale filtering procedure 
works for some values of A and /3; 



A 


/? 


a = -1/2 


a = 1/2 


P = A 


2 


1.6 


0<B< 4.57 


0<B< 1.52 


< S < 0.83 


1.5 


1.6 


< B < 1.35 


0< B < 0.54 


< B < 0.25 



As one can observe in the table, the reliability of the time-scale filtering procedure increases with 
respect to pattern 4 if we increase the influence of the decay-criterion. 

Pattern 5: For pattern 5 the situation is similar as with pattern 4. The edge at /? will - as 



0.2 0.4 0.6 



_/^-^ 



l! 0.4 Qi 0.8 



(a) a = 1/2 



(b) a = -1/2 



Figure 7: Numerical analysis of Qq for pattern 5. The figures display level curves of the surface 
Qo{A,B,l3) = for P(n,m) given by (1) with (a) a = 1/2 and (b) a = -1/2. 
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(a) a = 1/2 



(b) a = -1/2 



Figure 8: Numerical analysis of Qi for pattern 5. The figures display level curves of the surface 
Qi{A,B,l3) = for P(n,m) given by (1) with (a) a = 1/2 and (b) a = -1/2. 



before - 'push' the other maxima-lines away. However, since the 'long' maxima-line now approaches 
(and not 1 as for pattern 4, fig.3(h,j)), we expect that the distance-criterion will be strengthened. 
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Let as before Qi{A, B, (3) denote the minimum of the Ihs. of (4) for < ,82 < s*, $2 = 2si, j = 
and < A < A* . Given i?, /3 one can find A* by solving (2). Let Qq{A, B, /?) denote the minimum 
of the Ihs. of (5) for < si < s*, .$2 = 2si, j = 1 and < yl < ^*. 

Fig.7 and fig. 8 show some level-curves of the surfaces Qi[A, i?, /3) = 0, i = 0, 1 for P(n, m) de- 
fined in (1) with a = ±1/2. Given the pair [B, (3) in the i?/3-plane the point [B, (3, A) on the surface 
Qi{A,B, 13) = represents the maximal value of A such that the procedure successfully connects 
mod- max along the maxima-line it for i = 0, 1. 

One may observe that the reliability of the decision function P(-, •) increases if we increase the 
weighting of the distance-criterion. Indeed; Qi{A,B, (3) > 0, i = 0, 1 for 'almost' all 0<A<A* if 
P(n, m) = A(n, m) or if P(n, m) is given by (1) with a = 1/2. On the other hand, if a = —1/2 one 
can roughly say that one must have ^4 < 1 to ensure that both Qi > and Qq > 0. 

Pattern 6: Fig. 9 and fig. 10 show a few examples of the analysis of the decision- function with 
respect to pattern 6. Similar to the previous patterns let Qi{A, B , j3) denote the minimum of the 
Ihs. of (4) wrt. S2 = 2si, < si < s* , j = and A> A*. Given {B, /3) one can find A* by solving (2) 
(If A = A* the maxima-line at coarse scales trifurcates). Let Qq(A,B, (3) denote the minimum of 
the Ihs. of (5) wrt. S2 = 2si , < S2 < s* , j = 1 and A>A*. 

Fig. 9 and fig. 10 show some results of the analysis of pattern 6 with B = kA for 0<k<l and 
0.1 < /3 < 0.9. If /? is close to either or 1 one can in practice model the signal as pattern 1, pattern 
2 or pattern 3, depending on mutual relations between k and /?. For large values of (3 and k>l 
the maxima-lines become unstable making numerical analysis difficult. 

The surfaces Qi{A, B, f3)=0, i=0, 1 represent the minimal ^o>0 such that Qi{A, B, (3)>0 for all 
^>0. Fig. 9 shows some level-curves of the surface Qi{A,B, I3)=0 for P(n,m) given by (1) with 
a= lb 1/2. Fig. 10 shows the some level-curves of the surface Qq{A,B, f3)=0 for P{n,m) given by 
(1) with a = - 1/2. Numerical studies suggest that Qo{A, 5, /3) > for all 0.1 < /? < 0.9, 0<B <A 
and ^* < ^ if a = 1/2 or P = A, except possibly for A^A*. 

One may observe in fig. 9 and fig. 10 that increasing the influence of the distance-criterion in- 
creases the reliability of the scale-space filtering procedure for pattern 6. This is caused by the edge 
at < /3 < 1 which 'pushes' the other maxima-lines away in different directions and hence strength- 
ens the distance-criterion. In addition, the edge at /? has negative intensity B which disrupts the 
decay-coefficient of the wavelet transform, and hence the accuracy of the decay-criterion. 
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(a) a = 1/2 



(b) a = -1/2 



Figure 9: Numerical analysis of Qi for pattern 6. The figures display level curves of the surface 
Qi{A,B,l3) = for P(n,m) given by (1) with (a) a = 1/2 and (b) a = -1/2. 
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(a) Q = -1/2 

Figure 10: Numerical analysis of Qq for pattern 6. The figure display level curves of the surface 
Qo{A,B,l3) = for P(n,m) given by (1) with a = -1/2. 



3.4 Summary and discussion of 1-D time-scale filtering procedure. 



We have studied the suggested time-scale filtering procedure from Sect. 3.1 with respect to the six 
model-patterns in Sect. 3. 2. 

Connecting mod-max across time-scale plane for the first two patterns is trivial. For the other 
patterns the overall accuracy of the decision-function P(-, •) in general increases if we use a suitable 
weighting between the distance- and decay-criterion. For pattern 3 and pattern 4 the accuracy of 
the decision-function P(-, •) benefits from increasing the influence of the decay-criterion. For the last 
two patterns we saw that increasing the influence of the distance-criterion increases the accuracy of 
the time-scale filtering procedure. It is clear that one should take specific features of the concrete 
signals into account when defining the mutual weighting between distance- and decay-criterion. 
When analyzing the medical ultrasound images in Sect. 5. 3 we use a=0. 

In Sect. 5.1 we will evaluate the performance of the proposed time-scale filtering procedure with 
respect to medical ultrasound signals. 

The results from the previous sections are useful also for more complicated signals. For instance 
the signal f (t) = p(t)+Ap{t—l)+B' p(t— 7)— A' p{t—9)+p{t— 10) is locally similar to pattern 3 near 
the point and pattern 4 near the point 10. Due to fast decay of the wavelet ipsit) the maxima-lines 
of the 'local' patterns in f{t) and their exact equivalent are practically identical for < s < 3. It 
makes little sense to analyze the signal at scales s > 3 for the purpose of edge-detection. Our analy- 
sis of actual medical ultrasound signals indicates that such signals often are composed of (distinct) 
patterns in this manner as illustrated in fig. 2. 

Real medical ultrasound signals are of course 'smooth'. The corresponding edges can be mod- 
eled by smoothed versions of the 'ideal' patterns in Sect. 3. 2. If / is any of the ideal patterns say, 
then we can model the smoothed patterns as h = f * g^ where ga{t) = ,} e~* '^"^ ,0" >0, is the 
Gaussian smoothing kernel. 

With some natural modifications the suggested time-scale filtering procedure applies to the 
smoothed patterns as well. We need to use denser (i.e. non-dyadic) scales; more precisely we 
need to choose S2,si such that yjs^^+a^ = 2^Js\ + a^. If we in addition change ••• — 1/2 to 



1/2 






j^^ in the expression for the decay-criterion D(n,m), the results from Sect. 3. 3 apply 
directly (i.e. with the same A,B,/3). 
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3.5 Decision-function and space-scale filtering procedure for 2-D signals. 

We turn the attention towards a space-scale filtering procedure for connecting mod-max across the 
space-scale plane of 2-D signals. In images mod-max are chained together to form boundary- curves 
corresponding to boundaries of objects, and the natural space-scale plane object is the maxima- 
surface. The behavior of a maxima-surface is more subtle than for a maxima-line. For instance one 
can not assume that a maxima-surface does not branch. The target with the space-scale filtering 
procedure is to decide which mod-max propagate toward each of the boundary-curves at the finest 
scale. 

The suggested decision-function for 2-D signals is - with some natural modifications - simi- 
lar to the 1-D decision-function defined in Sect. 3.1. Let S2 and si be two scales S2>si, and let 
{n, S2)={{ni,n2), S2) and (tti, si)=((mi,m2), si) be two mod-max. We define the decision-function 
as; 

P(n, m) = A(n, m)D{n, ?TT-)Angle(n, m) 

where; 



A(n, m) — cxp( — |n — m\s-^ " 
D(n, m) ^ cxp 



^^ J^/(n,.2)| ^^-i^2 



\Wfim,si)\ si 

Angle(n, m) = cxp(— |^/(n, S2) — Af{m, si)\). 

Af{u, s) is the direction of the wavelet transform at the mod-max (u, s), and is given by 

f arctan ^©^ Wr{u, s) > 

Af(u s) = < wp{u,s) J \ , J - 

I 7r + arctan ^^^j^ VF/^(u,s)<0. 

The essence of the space-scale filtering procedure is - as in 1-D - to connect mod-max which 
maximize P(n, m). 

In particular, if one consider 2-D patterns obtained by a tensor product of the 1-D patterns, 
then the analysis in Sect. 3. 3 extends to the 2-D procedure. 

Given two scales S2 > si the procedure requires that for each n G Mf{s2) one has to compute 
P{n,m) for every m G Aif{si). One can reduce the computational efforts by computing P(n, m) 
only for m G Aif{si) within a given distance (4si should suffice if we use dyadic scales) of n. 

In the next section we will see that for our purpose it is sufficient to compute P{n, m) for only 
a few mod-max at scale S2 as well. By using these simplification one will effectively reduce the 
computational complexity of the space-scale filtering procedure for images. 

4 Edge-detection. 

The uncertainty-principle states that in a noisy signal there exists no single scale such that one can 
detect and simultaneously localize only the significant edges, [3]. However, by using the maxima- 
lines one can turn the uncertainty-principle into our advantage. At coarse scale one can detect the 
significant edges in the signal, and at fine scales one can localize them. Since maxima-lines connect 
coarse- and fine-scale information one can effectively improve edge-detection in noisy signals. 

In this section we discuss a multi-scale wavelet-based edge-detection scheme. The method is 
based on the suggested time-scale filtering procedure and a perceptual criteria similar to that in 
[6]. This approach has shown promising results for detecting and localizing (low-contrast) edges in 
medical images. 
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4.1 1-D edge-detection. 

Assume that the wavelet transform is computed at scales Sj for j = 1, ..., J, and let Mf{si) denote 
the mod-max at the finest scale si. For a mod-max a £ Mf{si) let ia denote the maxima-line 
which contains a. To each a assign the number; 

R{a)= Yl \Wf{n,s)\. 

If need be one can introduce a weight-function /i(s) in the expression for R{a), for instance to favor 
fine/coarse scale information. The final step of the edge-detector is to decide a threshold T > 
and find the mod-max a G A4f{si) such that R{a) > T. These mod-max are used to represent the 
significant edges in the signal. 

4.2 2-D edge-detection 

Assume that the wavelet transform is computed at scales Sj for j = 1, ..., J and let Cj denote the 
boundary-curves at scale si. 

To each maxima-curve Cj we assign the following quantity; 

S{cj) = (length of Cj) • av{{R{a) : a £ Cj}), (7) 

where av(-) denotes the average-function. The final step of the edge-detection is to identify the 
curves Cj such that S{cj) > T, where T>0 is a threshold. These boundary-curves is used to 
represent the boundaries of significant edges in the image. 

The quantity av({i?(a) : a £ Cj}) reduces the influence of a few false connections made by the 
space-scale filtering procedure along a boundary-curve Cj. Such false connections can occur for 
some 'critical' patterns (e.g. A ~ 1 in pattern 3). 

As mentioned in Sect. 3. 5 on can save computational efforts by reducing the number of rri's for 
which one has to compute P{n,m). Further reduction can be achieved by computing the average 
in (7) from only a 'few' randomly selected maxima-lines from each maxima-curve Cj. Experiments 
indicate that this simplification has little influence on the values of S{cj). 

5 Experiments and results. 

We test the edge-detector on both phantom and actual medical ultrasound signals. By studying the 
algorithm on phantom-signals we compare the Figure of Merit [19] with the Canny edge-detector. 
The visual performance of the edge-detector is analyzed with respect to medical ultrasound signals. 
This analysis focuses on some well-known problem features of medical ultrasound images; such as 
low-contrast edges and speckle-noise. 

Before studying the edge-detection algorithm we analyze the reliability of the suggested time- 
scale filtering procedure on actual signals. The main question is whether the results from Sect. 3 
also apply to actual medical ultrasound signals. 

5.1 Accuracy of time-scale filtering procedure on medical ultrasound signals. 

To quantitatively evaluate the performance of the suggested time-scale filtering procedure from 
Sect. 3.1 we compare the output of the procedure with the exact maxima-lines of the wavelet trans- 
form. Assuming that the procedure connects the mod-max (n, S2) to (m, si) and n=la{s2) say, then 
we will study if m=ia{si)- The test-signals used to evaluate the procedure are 1-D cross-sections 
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of medical ultrasound images. 

In the analysis we study the probability of a false connection, and between which scales an 
error is most likely to occur. For the purpose of edge-detection the severity of a false connection 
increases at coarse scales. We also study the average displacement (in pixels) of a mod-max caused 
by a false connection. If this number is small it indicates that an error has little influence on the 
visual output of the edge-detector. 

In the experiments we use scales s = 2-', 1 < j < 5. The control-parameter is chosen as 
a = { — 1/2,0, 1/2}^. To obtain the actual maxima-lines we use the edge-focusing algorithm [4]. 

The results of the quantitative analysis of the time-scale filtering procedure is displayed in 
Tab. 1-3. Each column represents the average performance of the procedure applied to every hor- 
izontal cross-section in the corresponding image. The images are displayed in fig.l2(a,b,c) and 
fig.llb. 

The analysis of pattern 6 in Sect. 3. 3 showed that the accuracy of the suggested time-scale 



Signal: 


Tumor 


Blood - 

vessel 


Liver 


phantom 


False connections (%): 


4.8 


4.7 


6.2 


4.2 


Spatial error in pixels: 


12.4 (2.5) 


12.9 (2.5) 


12.2 (2.8) 


20.3 (7.4) 


false connections between 
scale 32 and 16: 


2.5 


4.4 


8.5 


4.2 


false connections between 
scale 16 and 8: 


5.4 


5.2 


9.0 


8.0 


I'alse connections between 
scale 8 and 4: 


6.2 


5.7 


8.2 


5.0 


I'alse connections between 
scale 4 and 2: 


4.8 


4.5 


6.0 


1.7 



Table 1: Analysis of the time-scale filtering procedure in Sect. 3.1 on medical ultrasound signals. 
The control-parameter between distance- and decay-criteria is a=l/2. The number is parenthesis 
on the 2. row is the spatial error in percent relative to size of the signal. 



Signal: 


Tumor 


Blood - 

vessel 


Liver 


phantom 


False connections (%): 


4.0 


4.2 


4.8 


3.8 


Spatial error in pixels: 


11.5 (2.2) 


13.2 (2.6) 


12.6 (2.9) 


19.2 (7.0) 


I'alse connections between 
scale 32 and 16: 


1.9 


4.3 


8.6 


3.3 


I'alse connections between 
scale 16 and 8: 


4.3 


5.3 


7.9 


8.0 


I'alse connections between 
scale 8 and 4: 


5.2 


4.8 


5.4 


4.1 


I'alse connections between 
scale 4 and 2: 


4.1 


3.9 


4.7 


1.6 



Table 2: Analysis of the time-scale filtering procedure in Sect. 3.1 with control-parameter a=0. 



filtering procedure increases as a increases from —1/2 to 1/2. Speckle noise can often locally be 
modeled as pattern 6 (with A^B^l and /3«2)- Hence we expect to observe a similar increase in 
accuracy for actual medical ultrasound signals. In Tab. 1-3 one can indeed observe such increase 



^Note that in the analysis of the patterns in Sect. 3 we had si < 1 in the expression for P(n, m) in (1). In particular 
Since we are using scales s — 2-' , 1 < j < 5 to analyze actual medical ultrasound signals we 



„l/2 



<S? 



< S 



-1/2 



will necessarily have Sj^ 
distance-criteria. I.e. a 
a=l/2. 



1/2 



> sY > s- 



-1/2 



This causes that a large a favors decay-criteria and a small favors the 
1/2 in the preceding section 'corresponds' to a = 1/2 in this section, and visa verse for 
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Signal: 


Tumor 


Blood - 

vessel 


Liver 


phantom 


False connections (%): 


3.9 


4.0 


4.6 


3.7 


Spatial error in pixels: 


11.7 (2.3) 


13.5 (2.6) 


13.3 (3.1) 


19.3 (7.0) 


I'klse connections between 
scale 32 and 16: 


1.9 


4.3 


8.6 


3.3 


li'alse connections between 
scale 16 and 8: 


4.4 


5.3 


7.8 


7.8 


I'alse connections between 
scale 8 and 4: 


5.0 


4.5 


5.2 


4.0 


I'alse connections between 
scale 4 and 2: 


4.0 


3.7 


4.5 


1.6 



Table 3: Analysis of the time-scale filtering procedure in Sect. 3.1 with control-parameter a= — 1/2. 

of the accuracy of the procedure on medical ultrasound signals (recall that a=l/2 in this analysis 
corresponds to a= — 1/2 in the analysis in Sect. 3, etc). A closer inspection indicates that false 
connections in between maxima-lines corresponding to speckle-noise is the typical problem for the 
time-scale filtering procedure on ultrasound signal. Since speckle-noise contains information which 
is unwanted from an edge-detection viewpoint, such false connections will only be a minor problem. 



5.2 Results of edge-detection on phantom ultrasound images. 

To study the performance of the edge-detection algorithm suggested in Sect. 4. 2 we compare it with 
the well-known Canny edge-detector. We believe that comparing the Canny edge-detector with ours 
demonstrates some differences between multi-scale and single-scale wavelet-based edge-detection. 

For the Canny edge-detector with hysteria we have used thresholds Tiow = 0.1-max |VF/((-, •), s)| 
and Thigh = 0.3 • max |VF/((-, •), s)|. The scale is chosen equal to the finest scale used by our 
algorithm. 

To quantitatively analyze the edge-detectors we study the phantom medical ultrasound signals 
shown in fig. 11. The image in fig. 11(a) shows the noise-free version, and fig.ll(b,c) shows the 




(a) 



(b) 



(c) 



Figure 11: Phantom medical ultrasound images, (a) Noise-free image, (b,c) images with speckle- 
noise. 



images with various speckle-noise. The phantom images are generated by using a region image as 
a variance field for the diffuse scattering. The images are made by convolving the ideal image with 
a parametric point-spread function and adding white Gaussian noise, and are discussed in details 
in [20, 21]2. 



The images are kindly provided by SINTEF MedTech. The author thanks T. Lang0 for his assistance. 
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Signal: 


proposed edge-detector 


Canny 


phantom 1: 


0.56 


0.53 


phantom 2: 


0.54 


0.47 



Table 4: The Figure of Merit of the edge-detection algorithm in Sect. 4. 2 and the Canny edge- 
detector at scale 4. 



Signal: 


proposed edge-detector 


Canny 


phantom 1: 


0.54 


0.53 


phantom 2: 


0.53 


0.51 



Table 5: The Figure of Merit of the edge-detection algorithm in Sect. 4. 2 and the Canny edge- 
detector at scale 6. 



We use the well-known Figure of Merit (F0M)[19] to quantitatively evaluate the edge-detection 
algorithms. The FOM is the number defined as; 



FOM 



1 



maxjrip, n^} 






where Up and rid are the number of resp. true and detected edges, and di is the minimal distance 
from a detected to a true edge. 7 is a constant used to penalize displaced mod-max. To be consistent 
with others we use 7=0.11 [8, 22]. The FOM should ideally be close to 1. 

The values of the Figure of Merit of the two edge-detectors are displayed in tab. 4- 5. 

5.3 Result of edge-detection on real ultrasound images. 

In this section we discuss the results of the suggested edge-detection algorithm on real signals. For 
this purpose we study the edges detected by the algorithm in several actual medical ultrasound 
images. The goal is to study if the algorithm detects the significant edges in the images without 
including irrelevant information such as noise. In addition we discuss how one can obtain a threshold 
as well as the computational complexity of the algorithm. 

Fig. 12 shows the ultrasound images which are discussed in this section. Each image contains 




(a) (b) (c) 

Figure 12: (a) : Tumor-image, (b) : Liver-image, (c) : Blood-vessel image, 
some typical features of medical ultrasound images which cause difficulties for edge-detection. 
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Fig. 12(a), Brain-tumor image: The interesting edges are the outer boundary of the tumor and 
the cortex (white 'horizontal' object to the south-west of the tumor). The voids within the tumor 
are mainly interesting in order to study how the edge-detector represents fine details. 

Fig. 12(b), Liver image: The interesting edge is the boundary of the liver. The boundary is the 
horizontal line in the upper third in the image. This image interesting to study because of the low 
CNR of the boundary. 

Fig. 12(c), Blood-vessel image: In this image we want to detect the boundary of the blood-vessel. 
Some parts of the boundary of the blood-vessel have low contrast making them difficult to detect. 
This problem is typical for ultrasound images of blood-vessels. 

Edge- detection and noise-removal: Fig.l3(a,b,c) shows some results of our edge-detection 




vj. 




(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 13: Results of the edge-detection algorithm of the images in fig. 12. In figure (a) - (c) we 
display the results of the edge-detection algorithm studied in Sect. 4 . In figure (d)-(f) we display 
the results obtained with a Canny edge-detector. (a,d) : Tumor image, (b,e) : Liver image, (c,f) : 
Blood-vessel image. 



algorithm applied to the images in fig. 12. Most of the significant edges which we would like to 
detect are present in the representations. One may also observe, in particular in fig.l3(a,c), that 
the representations contain few edges which can be ascribed to noise. The representations in 
fig. 13 show promising results for detecting low contrast edges and fine detail edges. However, the 
boundary-curves of such edges are vulnerable to fragmenting (which make them difficult to detect). 
This phenomenon occurs for instance underneath the tumor (low contrast) and to some degree 
inside the tumor (fine details). The fragmenting - in particular for low-contrast edges - appears to 
increase due to speckle-noise. One can reduce this problem by connecting fragmented boundary- 
curves (e.g. by using coarse scale information as suggested in [6]). Preliminary studies have shown 
that using this as a complimentary tool improves edge-detection of low-contrast and fine detail 
edges. We use no synthesizing of the boundary-curves in this study. 
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Fig.l3(d,e,f) shows the result obtained by using the Canny edge-detector (with hysteria). The 
scales is chosen equal the finest scale used by our algorithm. One can observe that the edge-detector 
scheme studied in Sect. 4. 2 in general yields better representations of the significant edges in the 
images, and includes less noise. 

Threshold: The edge-detector assigns to each boundary-curve Cj a value S{cj) representing both 
its length and value of wavelet transform in both time- and scale-space. Values of S corresponding 
to noise will therefore be concentrated at low values, while those corresponding to edges will be 
distributed discretely at high values. The result is that there is often a separation between values 
of S corresponding to noise and edges. A similar phenomenon was observed in [6]. In addition, 
since values of S corresponding to significant edges are often distributed discretely it is possible to 
find representations which are perceptually optimal. 

Computational complexity: The computational complexity of the edge-filtering procedure in 
Sect. 4. 2 is negligible. In what follows we focus on the computational complexity of computing the 
wavelet transform and the space-scale filtering procedure from Sect. 3. 5. 

To reduce the computational complexity of the space-scale filtering procedure we use the tricks 
discussed in Sect. 3. 5 and Sect. 4. 2. In other words we find the maxima-lines of only a few randomly 
selected mod- max along each boundary-curve, and use a sliding window (of size 4 -scale x 4 -scale) to 
decide candidate connections at the next finer scale. We use the fast Fourier-transform to compute 
the wavelet transform. 

The space-scale filtering procedure has been implemented in MATLAB version 7.9.0 (R2009b) 
on a laptop with a Pentium 3 (2 GHz) processor. Tested on a 362 x 512 image the total processing 
time is 4.6 seconds if we find the maxima-lines of 50% of the mod-max. If we find only 10% 
randomly selected maxima-lines the total processing time is 1.9 sec. 

Computationally both computing the wavelet transform at different scales and connecting mod- 
max across scales is well-suited for parallel-computing. It is expected that this should allow us to 
further reduce the total processing-time of the scale-space filtering procedure. 



6 Conclusion 

We suggest a fast multi-scale edge-detection algorithm for medical ultrasound signals. The edge- 
detector is based on the continuous wavelet transform, and takes advantage of the behavior of 
maxima-lines in the space-scale plane. Our algorithm shows promising results for detecting the 
significant edges in medical ultrasound signals. 

The key-point of the edge-detection scheme is a fast space-scale filtering procedure for medical 
ultrasound signals. The procedure is designed to accommodate a set of patterns modeling intensity- 
changes which are typical in such signals. This allow us to manage and reduce the ambiguity intro- 
duced when using the multi-scale wavelet transform to detect edges in a medical ultrasound signal. 
The space-scale filtering procedure connects mod-max which belong to the same maxima-line by 
computing the wavelet transform at only a sparse set of scales. For our medical ultrasound signals 
we typically need to compute the wavelet transform at the four scales s = {4,8,16,32}. In the 
popular Edge-focusing algorithm [4] one computes the wavelet transform at a linear set of scales 
(As = 1/2) to trace mod-max in the space-scale plane. Our proposed space-scale filtering proce- 
dure will therefore considerably reduce the computational efforts needed to decide which mod-max 
belong to the same maxima-line. 

Based on the space-scale filtering procedure we study a novel edge-detection algorithm for med- 
ical ultrasound signals. The edge-detector assigns to each candidate edge a value representing the 
behavior of its corresponding maxima-line in the space-scale plane. Experimental results indicate 
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that the edge-detector finds both high- and low-contrast edges in medical ultrasound images and 
yet disregards noise. Compared to the Canny edge-detector we observe that our algorithm finds 
the significant edges, and includes less redundant information such as noise. Our study has showed 
that one can construct multi-scale edge-detection algorithm for medical ultrasound signals which 
meets the demand of speed, but not on the expense of reliability. 

We also expect that using multi-scale information opens new possibilities for further improve- 
ments of edge-detection in medical ultrasound images. One possibility is to dynamically choose the 
finest scale based on local automatic scale selection as described in e.g. [16]. 
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